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laid near the turn of the century by the discovery of Fe 
and Ni carbonyls, culminating in the potthuious publication 
of a body of data on the synthesis, stability and reactions 
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A promising Method (or extracting oxygen fro* ileenite 
is through the hydrogen-reduction procee* (f ig. 4 ) . the 
relevant reaction Is: 
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THEORETICAL PREDICTIONS OF VOLATILE BEARING PHASES AND VOLATILE RESOURCES 
IN SOME CARBONACEOUS CHONDRITES 


Jibamitra Ganguly 
Department: ot Geosciences 
University of Arizona, Tucson, AZ 3 57 2 1 

- u r e no r a * , s a x q na 
Department of Geology 
Brooklyn Cc liege, Brooklyn, *.Y 11280 


Abstract 

Darbonaceous cr.cnarites are usually 
believed to be the 1 ' '"primary constituents of 
near-Earth asteroids ana rnobos tia 
Dierr.os, and are potential resources of 
fuels wnich may be exploited for future 
planetary missions. In this worK we have 
calculated the nature and abundances of 
the major volatile bearing and other 
phases, including the vapor phase, feat 
snould form in Cl ana C2 type carbonaceous 
chondrites as functions of pressure ana 
temperature. The results suggest that 
talc, ancigorite x magnesite are the major 
volatile bearing phases and are stable 
below 400°C at 1 bar in these cnoncntic 
compositions. Simulated heating of a 
kilogram of C 2 enonente at ;ixec bulk 
composition between 400 and 800 v r at 1 bar 
vieids aoout 135 gm of volatile, wnich is 
naoe primarily of HjO, H 2 O, CH 4 , CO-i and 
CO. The relative abundances of these 
volatile species change as functions of 
temoe r a ture, and on a molar basis, Ho 
becomes the most dominant species above 
550°C. In contrast, Cl chondrites yield 
about 306 gm ct volatile under the same 
condition, which consist almost completely 
of 60 wt % H and 40 vt% CO 2 * Preliminary 
kinetic considerations suggest tnat equi- 
librium dehydration of hydrous phyllosili- 
cates should be attainable within a few 
hours at 600°C. These results provide the 
framework for further analyses of the 
volatile and economic resource potentials 
of carbonaceous chonantes. 

In t roduct ion 

The near-Eartr, asteroids are poten- 
tial resources of volatiles which can be 
used as propellants and life supporting 
purposes for planetary missions. As 
emphasized by Lewis and Lewis (l), these 
asteroids are at times. the nearest bodies, 
to Earth, and many of them can pass 
between Earth and Moon. These authors 
have also shown that within any arbitrary 
three-year period, there are roughly 300 
and 90,000 launch opportunities to respec- 
tively kilometer and 100 -meter sized 

asteroidal bodies as compared to 39 lunar 
launen windows. These tacts and the 
potential volatile resources make the 
asteroids important candidates tor 
detailed analyses tor supporting future 
planetary missions. 

Cost-effective or energy-et t ic len t 
extraction of volatiles from the asteroids 
require knowledge of the nature ana modal 
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-ouncances of tne constituent minerals in 
wnicn tr.e volatile •: cm pc cents are struc- 
turally oounc. There nas neon several 
studies on reflectance cpectroscopv of 
near-Eartr. asteroids to characterize the 
nature of the volatile '-earing minerals 
: 2 , 3 ) . v.hi le these studies report absorp- 
tion band near 3 -m suggesting presence Of 
hydrous phy lios 1 1 icates , -ore ^recise 
characterizations were ■qyi vocal. An 

obvious alternative approach to the reso- 
lution :r the problem is to examine 
directly tne asteroidal materials. 

The spectra i and density cnaracter- 
istics strongly suggest that 2 significant 
fraction of the near-Eartr. ssteroias are 
^ a d e o i carbonaceous j r. 0 n a r 1 1 e s - 1,4). 
The most volatile non of these meteorites 
are wnat are common iy known 2 s Cl or (Cl) 
and 0 2 (or DM) -lasses. However, the 
;rai n size of the volatile no a ring phases 
are otten too small l IC>0 t: I >00 A) (5) to 
permit identification unaer optical micro- 
scope, ana are best observed jnaer high 
resolution transmission T.-;^ctron micro- 
scopy * H STEM ) Dolensky vo ^cEween (5) 
have recently presented a summary ot vola- 
tile bearing phases identified m the 
carbonaco us cnona rites. 

The primary objective c f this work is 
to carry cut theoretical : a i cu 1 a t ions to 
predict modal abundances int compositions 
of tne na jor mineral pnasos, along vith 
the abundance and composition the co- 

existing vapor phase, that :ouid develop 
tn the bulk compositions :■ t D 1 and C 2 
choncrites as functions ot pressure ,n<i 
temperature. The results wo-jio L rovine 
the framework ror engineer mo or signs tor 
..tne extraction of volatile coipononts fr^m 
_asteroids as well as Pnobos ana Deimos, 
_ t he two small natural satellites of Mars, 
which are also likely to be mane of car- 
bonaceous chondrites (I). 

Theoretical Method 

Principles 

According ro Du hen's theorem in 
--.classical thermoriynam ics the equi- 

librium state of a closed system ii.e., 
system or fixed composition and mass) is 
completely determined ' i f my two variables 
are fixed regardless ot whether these are 
intensive, extensive, or a combination of 
both. Ganguly ana Saxena < 7 J have recent- 
ly reviewed tne various methods by which 
one may carry out the actual computation 
of the equilibrium assemblages. The 
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Table 1 Bulk chemical 

compositions of Cl and C2 
carbonaceous chondrites 



Cl 

C 2 

Si 

10.40 

12.96 

Ti 

0.04 

0.06 

Al 

0.84 

1 .17 

Cr 

0.23 

0. 29 

Fe 

18.67 

21.56 

Mn 

0.17 

0 . 16 

Mg 

9.60 

11.72 

Ca 

1.01 

1.32 

Na 

0.55 

0.42 

K 

0.05 

0.06 

P 

0.14 

0.13 

Ni 

1.03 

1 . 25 

Co 

0.05 

0.06 

S 

5.92 

3.38 

H 

2.08 

1.42 

C 

3.61 

2.30 

0 

-45.61 

41.74 


method used in this study obtains the com- 
positions and abundances of the equilib- 
rium phase assemblages by minimizing 'the 
Gibbs Free Energy (G) of the system at 
fixed P-T conditions. The minimization is 
constrained to conserve the various ele- 
mental masses of the system through the 
method of Lagrangian multipliers < 7 ) . 
Because of the extremely small grain size 
of the phy 1 los i 1 ica t es ODserved in carbo- 
naceous chondrites, one should consider 
the effects of surface free energy on the 
formation of these phases. As a first 
approximation, we have, however, ignoced 
this effect. Therefore, the actual equi- 
librium temperatures of volatile bearing 
phases could be lower than those predicted 
by our calculations. 

System, Phases and Data Base 

Table 1 shows the bulk compositions 
of Cl and C2 carbonaceous chondrites. 
These compositions are based on analyses 
of a number of fragments ot Crgueil and 


Table 2 Phases in the system 
Mg-Fe-S i-C-H-O-S considered in this work 

The sources of t hermochem lea 1 data are shown by reference 
numbers within parentheses. The Heat ot Formation of the Fo-enn 


members of the 

hydrous phy 1 losi 1 icates are derived in 

this work. 


Anhydrous Silicates and Oxides 


01 i vine 

{ Mg , Fe ) iu 4 

( 3 ) 

Orthopyroxene 

(Mg,Fe)Si0 3 

( 3 ) 

Pe ri clase 

( Mg , Fe JO 

(23) 

Ouart z 

Si °2 

(17) 

Cristobal ite 

Si0 2 

(17) 

Hematite 

Fe 2 ° 3 

( 23 ) 

Magnet ite 

Fe 3 0 4 

< 23 ) 

Hydrous and Carbonate Phases 


Anthophy 11 ite 

( Mg , F e ) 7 S ig 0 2 2 i OH 1 ~> 

(8,27) 

Talc 

(Mg,Fe) 3 Si 4 O I 0 (OH) 

( 3 ) 

An t igor i te 

(Mg,Fe) 4 Q3i34Og5(0H)^2 

( 8 ) 

Chrysot i le 

(Mg,Fe) 3 Si 2 0 5 {0H ) 4 

( 8 ) 

3ruc ite 

( Mg , Fe ) ( OH ) j 

(17,24) 

Magnesite 

(Mg,Fe)C0 3 

(17,25) 

Sulfides ana E 

lene n t s 


Tro l 1 i te 

FeS 

(17) 

Py ri te 

FeS 2 

( 17 ) 

T ron 

Fe 

1 26 ) 

Sulfur 

S 

(17) 

Graph ite 

C 

E 17 ) 

VaDor Phase 

C-O-H-S 

(15,16,17) 

Species: H 2 O, 

CO 2 , CH 4 , CO, 0->, H ->S , S ^ , SO^ , COS 
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Murchison meteorites, respectively, ana 
are summarized in Dodd s22). It can De 
easily seen from this table that the sub- 
system Mg-Fe-S l-L'-H-O-S (MKCHOS) consti- 
tutes almost 96i by weight of toe bulk 
compositions ot Cl and C2 chonarites. In 
oraer to somewhat simplify t ne computa- 
tional problem, we have confined the G 
minimization calculation to this sub- 
system, but evidently the phases whose 
compos i t ions form outside ot this sub- 
system cannot be major constituents of 
these chonarites. 

Table 2 shows the list of phases that 
we have considered as possible crystal- 
lizing phases within tne above subsystem 
along with the sources of the thermochem- 
ical data* The choice of these phases 
have been guided by the puolished reports 
of the mineralogy of carbonaceous chon- 
drites whose compositions lie within the 
system MFCHOS {e.g., 2olensky and McSween 
(5)), The selection of thermochemical 
data in the system Mg-Si-Q-H is based on 
the critical analysis of Chatterjee (8,11) 
tor consistency with both calorimetric and 
exper unenta 1 Iy determined phase equilib- 
rium measurements. These data are similar 
to those of Berman ( 9 i . However, unlike 
those in Berman, the data are also consis- 
tent with high P-T experimental data 
involving distribution of Fe and Mg among 
coexisting silicates. 

The heat capacities (C p ) and entro- 
pies (5°) of the phases, if^not available 
in any of the recent systematizations of 
self-consistent data set, are estimated by 
stoichiometric summation of the properties 
of oxides and structurally analogous com- 
pounds, as discussed in a number of publi- 
cations (13,14). 

There are no thermochemical data for 
the Fe-end members of the phy 1 los i i lea tes 
listed in Table 2. The enthalpies of 
these phases are estimated according to 
the known relative properties of Fe- and 
Mg -end members of biotite, which is also a 
phy 1 los i 1 icate . The data derived by 
Chatterjee (3) and Robie Hemingway 

(10) for the Mg- and be -end members 
(phlogopite and annite, respectively) 
yield an enthalpy difference of -359.566 
KJ per mole of divalent cation at 1 bar, 
298 K. We have used these values to esti- 
mate the enthalpies of formation from 
elements of Fe-end members of phyllosill- 
cates . at 1 bar, 298 K from the available 
data for the Mg-end members. The results 
are as follows {KJ/mol of Fe). Fe-talc: 
-1606.80; Fe-ant igorite : -1 1 27.53 and Fe- 
chrysotile: -1094.87. The corresponding 

values of the free energy of formation 
from elements (aG?) for Fe-talc and Fe- 
antigorite at 1 bar, 29^ K are -4477,18 
and -4784.0 KJ/mol of Fe* , respectively. 

The validity of the estimated thermo- 
cnemical prope rt les of the Fe-end members 
of the phy llosil icates may be tested as 
follows. The Fe-end members of talc and 


antigorite are unstable. Insteaa, the 
stable phases are closeiy analogous com- 
pounds, namely, minnesotite, which has a 
stoichiometry :-c F^gS 1 36 c 86 ( 0H ] 26 3S 
comparea to that of Fe^TJiicOgQtOHJjQ for 
Fe-talc, ana g reena i i te ' w i tri a stoichiom- 
etry -f F 64361 32 ° 80 < 0H } 64 compared to 
FG48* 1 ’ i i i ,,J 8 5 ( OH ) 52 Fe-ant igon te. 

Anovicz’et ai. i!2) nave recently decer- 
minea L rr e of minnesotite and g reena lite 
from ^r.ase equilibrium o a t a . Their 

results, ^;.G f (minnesotite) = -4577.73 KJ 

3nd loj **' j reena 1 i te ) - -4799.50 K J , 

differ in the ngnt direction from the 
corresponding values of Fe-talc and Ke- 
an t igorite, considering the relative 
stabilities of the phases. 

The stereochemical environments of 
divalent cations in biotite are quite 
similar to those m talc, both being 2:1 
(TOT) layer silicates, and partly similar 
to those in antigorite, wnich is an 1:1 
(TO) layer silicate (28). Thus, these 
phyllosil icates should have similar 
differences _ between the enthalpies of 
formation of the Fe- and .Mg-end members. 
Further, since on the oasis of observa- 
tional data, the predicted compositions of 
the mineral phases in carbonaceous chon- 
drites are expected to be Mg- r i ch , the 
effects of errors tn the Fe-end member 
thermochemical properties and solid solu- 
tion model will be relatively small. This 
can be easily understood by considering 
the form of G-X curve of a stable solution 
and evaluating the effects of errors in 
above properties near Mg-terminal segment. 

The .jJjermodynamic mixing properties 
of the Fe“ - and Mg-end memoer components 
of the anhydrous silicates are taken from 
the sel f-ccnsistent summary of Chatterjee 
(8), while those of the pnyllos i 1 icates 
are assumed to be ideal for the lack ot 
any data. The latter assumption is likely 
to be approximately valid at T 2_ 600° C, 
since at these conditions biotite solid 
solution is found to behave approximately 
ideally (7). At lower temperature, the 
phyl losi 1 icates may deviate from ideal 
solution behavior, but pernaps not too 
dissimilarly to significantly affect their 
relative stabilities. Magnesite is also 
assumed to be an ideal solution of Fe-Mg 
components. As emphasized by Ganguly and 
Saxena (7), the computed dehydration con- 
ditions are not very sensitive to errors 
in activity-composition relations since 
the devolatilization equilibria are char- 
acterized oy large enthalpy changes. The 
iron oxides and sulfides are considered to 
be stoicmometric phases since the solid 
solutions of Mg are around a few percent 
at the temperatures of interest. 

The fluid properties are taken from 
Saxena and Pei (15,16) and Robie et al. 
(17). The “corresponding state method" of 
estimating P-V-T relation, as discussed 
by these authors, is extended to include 
the sulfur species. For mixing data, we 
have assumed that H 2 S behave as H 2 0, S 2 as 


m 




Mol/Kg _ Mol/Kg 



Temperature (K) 


5.00 - 



Temperature (K) 


Fig. 1 Calculated modal abundances of minerals per kilogram of 
cock in (a) Cl and ( b) C2 caroonaceous ehondritic bulk compositions 
as functions of temperature at 2 Kb. The steps in the calculations 
are shown by symools. 


- t r and SO 2 ana COS as CO-i , ana used the 
equivalent binary interaction parameters 
from saxena and Fe i (15). 

Resu 1 1 s 

At 1 bar, the free energy minimiza- 
tion calculations do not / 1 e 1 <4 any vola- 
tile bearing phase at T >_ 400 "C, At the 
present state of the t hermochem ica 1 data, 


are reliable at lower temperatures. In- 
deed, we find results at lower tempera- 
tures which do not seem to follow the 
systematics of the results obtained at 
nigher temperatures. Thus, we have calcu- 
lated the equilibrium pnase assemolages at 
1 pressure of 2 Kb so that the volatile 
bearing phases appear at higher tempera- 
tures. This procedure at least permits 
recognition of the major volatile bearing 
phases. Some of these peases may be re- 


-e co not 


that cur ca i.eu la t ions 
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Fig. 2 Ja leu laced compositions of the silicate phases in ci ana 
C2 carbonaceous chondrites as tunctions o f temperature at 2 Kb. 


placed by more stable assemblages at lower 
temperature, but should reappear during 
■eating experiment. The results, which 
are illustrated in Figure 1, snow that the 
primary volatile bearing phases in the Cl 
and C2 bulk compositions are talc, anti- 
gorite ± magnesite. Talc is stable to 
around 650°C at 2 Kb, and dehydrates to 
orthopyroxene ( OPx ) at higher temperature 
according to 7c + 01 - 5 OPx *- H^O. Anti- 
go rite and magnesite crystallize at T < 
350°c, but magnesite is absent in C2 bulk 
composition at T 400 V C. r'eS is a stable 
pnase at all temperatures investigated in 
this work. 've , however, did not consider 
the formation ot sulf3te or pyrrhotite. 
Magnetite cecomes a stable phase at 
T 500°C. Elemental sulfur found :n Cl 
cnondrite is invariably associated with 
pyrrhotite [Fei_ x Si (3). This suggests 
decomposition o‘f troilite {FeS) to pyrrho- 
tite and sulfur. 

Figure 2 shows the compositions of 
the f er romag nes i an phases as function of 
temperature at 2 Kb for the Cl ana C2 bulk 
compositions. Calculations at 1 bar also 
yield similar compositions. The observed 
compositions of these phases in the car- 
bonaceous chondrites vary widely in Fe/Mg 
ratio. One way to achieve this diverse 
i-e/Mg ratio is by aqueous alteration which 
seems to be a characteristic of the car- 
bonaceous chonarites. Fe*' fractionates 
very strongly into an aqueous phase rela- 
tive to the silicates (18). Thus, varying 
degrees of equilibration with an aqueous 
phase could lead to a spectrum of Fe/Mg 
compositions of the f erromagnesian sili- 


cates. Aqueous alteration may not be the 
sole reason, but is likely to be an impor- 
tant contributing t actor in the develop- 
ment of variable Fe/Mg ratio of the ferro- 
magnesian silicates in carbonaceous chon- 
j r i tes . 

The abundances of volatile species in 
equilibrium with the solids per kilogram 
of total mass of C2 chondrites are illus- 
trated in Figure 3. Figures 3a and 3b 
illustrate the molar ana mass abundances, 
respectively, as functions of temperature 
at 1 bar. It should be noted that the 
total masg of the volatiles do not change 
above 400 C, wnich means that under equi- 
librium condition C2 chondntic material 
wil^ completely aevoiatize if heated above 
400 c. The total yield of the volatile is 
about 14% of the initial mass of C2 chon- 
drites. The change in the abundance of 
various species above 400°C is due to 
homogeneous reactions within the fluid 
phase. Above 550 C, hydrogen has the 
highest molar concentration in the vapor 
phase. Hydrogen is a very good propellant 
and reducing agent which may thus be pref- 
erentially extracted by heating the C2 
chondntic material to T > SSO^c. 

In contrast to the C 2 chondrites, the 
volgtiles given off by Cl material above 
400 C consist essentially of H 2 0 and CO 2 
at 1 bar pressure. Per kilogram of mate- 
rial, the total yield of volatile is about 
306 grams of which nearly 60 wt% is H 2 O. 

The above results are subject to 
considerable uncer ta inties owing to the 
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Temperature (C) 



Temperature (C) 


Fig 

the 


3 Calculated equilibrium (a) molar and (b) 
volatile species per kilogram of C2 chondritic 


mass abundances 
material . 


of 


approximate nature of the thermochemical 
properties of the phy 1 los 1 1 ica te solid 
solutions, and also since we have ignored 
the effects of surface free energy on the 
formation of extremely fine grained crys- 
tals. Further, carbonaceous chondrites of 
any given type do not have homogeneous 
compositions, but instead probably repre- 
sent agglomeration of materials formed at 
different conditions at different parts of 
the solar system. However, these calcula- 
tions are still useful in providing an 
iaea of the relative stabilities, as well 


as upper stabilities of the volatile bear- 
ing phases and their average abundance. 
It is interesting to note that recent 
reflectance spectroscopic studies of car- 
bonaceous chondrites seem to corroborate 
the predictions made above as to the 
nature of the most abundant hydrous phases 
(Gaffey, pers. comm.). 

Kinetics of Devolatilization 

The equilibrium calculations pre- 
sented above provide a framework for the 
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7Tglc — 3Anth 4- 4Qz + 4H20 



Temperature (C) 


Fig. 4 Dehydration of talc as functions of time, temper- 
ature and extent of overstepping above the equilibrium 
dehydration boundary at 1 Kb. Volume fraction is that of 
talc remaining after reaction. Experimental data points, 
shown by asterisks, are from Greenwood {19}* 



1QQ 200 000 400 500 600 700 800 900 

Temperature (C) 


Fig. 5 Equilibrium relations of talc, 
antigorite and brucite along with other 
minerals in the system MgO-S 
Reproduced from Evans and Guggenneim (21). 


evaluation of the volatile as well as some 
mineralogical resource potential of Cl and 
C2 carbonaceous chondrites. The actual 
extraction of volatiles as a function of 
temperature and heating rate, however, 
depends on the kinetics of the devolatili- 
zation process. For kinetic analysis, it 
is essential to identify the minerals in 
which the volatiles are most likely to be 
structurally bound. The above results 
suggest that talc, antigorite and magne- 
site are likely to be the major volatile 
bearing phases in carbonaceous chondrites. 

Except for pure Mg-talc, there are no 
kinetic data for the other volatile bear- 
ing phases. The dehydration kinetics of 
talc was determined by Greenwood (19). At 
a given temperature, the rate o| reaction 
depends on the rate constant, K / and the 
extent of departure from the equilibrium 
boundary. This is expressed in transition 
state theory |20) according to the rela- 
tion Fate = K { 1-exp ( nAG/RT ) J , where n is 
a constant, commonly assumed to be unity, 
and AG is the Gibbs free energy of the 
reaction at the P,T condition of interest. 
Figure 4, which is constructed from Green- 
wood's experimental data, illustrates the 
dehydration kinetics of talc as functions 
of time, temperature and extent of over- 
stepping (AT) above the equilibrium dehy- 
dration boundary at 1 Kb. It is evident 
from this figure that essentially complete 
dehydration of talc can be achieved within 
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6 hours at about 160°C above the equilib- 
rium boundary. 

The stability fields of the pure Mg- 
end members of various phyllosilica tes of 
interest in this study are summarized for 

the condition of P tot al 3 P H20 in Figure 
5, which is reproduced from Evans and 
Guggenheim (21). According to our re- 
sults, talc and other phy 1 los l 1 ica tes com- 
pletely dehydrate at 400 C at 1 bar in 
Cl and C2 bulk compositions. This is com- 
patible with the phase relations shown in 
Figure 5. Thus, based on the results 
shown in Figure 4, one would anticipate 
that essentially complete dehydration ot 
all p hy 1 1 o s i 1 i c a t e s in Cl and C2 cbondrit— 
ic material is likely to be achieved ig a 
matter of a few hours at T 2. 600 c. 
Detailed studies of the devolatilization 
kinetics of both phyl los 1 1 icates and car- 
bonates as functions of temperature, fluid 
composition, and grain size are currently 
under way to evaluate the practical limi- 
tations of the extraction of volatiles 
from carbonaceous chondritic materials. .... 
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